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Resonant annihilation of long-lived massive colored particles through hadronic 

collisions 
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Hypothetical long-lived massive colored particles (MCPs or Ys) would be confined in colorless 
exotic strongly interacting massive particles (SIMPs) at color confinement temperature of Tc ~ 
200 MeV. Two long-lived MCPs form a bound state {YY) at collisions of two SIMPs. We study 
sensitivities of MCP annihilation to decay properties of resonances {YY), and binding energies or 
energy levels of exotic SIMPs. The {YY) formation is assumed to dominantly proceed through 
resonances of {YY) in this paper. We make a toy model of the effective cross section for YY 
annihilation. Abundances of SIMPs are then calculated for different sets of parameters specifying 
properties of {YY) resonances, binding energies of SIMPs, the initial abundance and the mass of 
MCP. Calculated relic abundances for respective SIMP species are 2 x 10~*-3 x 10~* times that 
of baryon. They can be much higher but cannot be much smaller than the previous estimate. 
The abundances can be consistent depending on parameters with the possible scenario that SIMPs 
bind to nuclei and subsequent exotic nuclear reactions reduce the primordial abundance of ^Li or 
enhance those of ^Be and/or B in the early Universe. A unique information on the quark-hadron 
phase transition in the early Universe may become available in future by elaborated studies on the 
annihilation process with light element abundances as observables. 

PACS numbers: 98.80.Cq, 14.80. Pq, 13.75.-n, 98.80.Es 



I. INTRODUCTION 

In some particle models beyond the standard model 
(SM). long-lived massive colored particles (MCPs) can 
exist. The models include a long-lived gluino in split su- 
persymmetry [J, 0] , the next-to-lightest supersymmetric 
particle (NLSP) gluino [HI, and NLSP squark [| both 
in a weak scale supersymmetry. In addition, extended 
theories with other kinds of colored particles have been 
considered 0, Such models will be tested at weak 
scale in the Large Hadron Collider (LHC). 

If long-lived MCPs exist at the temperature of color de- 
confinement, i.e., Tc ^ 200 MeV, they would be confined 
inside exotic heavy hadrons, so-called strongly interact- 
ing massive particles (SIMPs, or X particles) 

Dover, Gaisser & Steigman [l^ have originally esti- 
mated the relic abundance of SIMPs, and studied the 
possibility of exotic stable hadrons with limits from ter- 
restrial superheavy element searches. They adopted the 
thermal rate of MCP annihilation in collisions of SIMPs 
which is assumed to be the same as (or less than) that 
for nucleons, and av ~ 30 mb. The final abundance thus 
estimated was about 10~^^ times the baryon abundance. 

Nardi & Roulet have studied the possibility of new 
stable exotic quarks predicted by the Eg model. They 
concluded that the possibility is ruled out unless there 
is some mechanism for decay of such quarks. Also, they 
pointed out that the annihilation cross section used in 



Ref. [l^l is significantly overestimated. The first reason is 
that the cross section includes extra final states irrelevant 
to the MCP annihilation. The second is the de Broglie 
wave length of a two-MCP system is the maximum radius 
scale for the total annihilation cross section. It is never 
larger than the interaction radius of -^1 fm until the tem- 
perature decreases down to ^10 MeV (m/1 TeV)~^ with 
m the mass (to ^ 1 GeV) of the MCPs. However, the 
previous work [lol | used a cross section larger than 1 fm~^ 
above the temperature. 

Thermal relic abundances after the freeze-out of an- 
nihilations in the early Universe have been studied with 
various assumptions on annihilation cross sections [l3j . 
Theoretical estimates, however, predicted various values 
ranging over more than several order of magnitude de- 
pending on adopted cross sections [13 1. 

Kang, Luty and Nasri [3l studied aspects of annihila- 
tion of SIMPs at T - Tc. They showed that the annihila- 
tion proceeds through the following steps: the formation 
of bound states of two MCPs with large orbital angu- 
lar momentums, their transitions to more bound lower- 
energy levels, and the annihilation of MCPs. They sug- 
gested that the bound state formation can realize from 
multiple partial waves of initial relative angular momenta 
in the two-SIMP system. They then suppose a geometri- 
cal cross section for annihilation at temperatures T ~ Tn . 

The final abundance of SIMP has been estimated 
to be 
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where Nx is the number density of SIMP, s ~ 
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27r^g*s7^^/45 is the entropy density with g^,s ~ 10 the 
total number of effective massless degrees of freedom in 
terms of entropy [l^ below the phase transition of quan- 
tum chromodynamics (QCD), is the total number of 
effective massless degrees of freedom in number ^Ej, a 
is the annihilation cross section of the SIMP, Tb is the 
temperature of the Universe at which the colorless SIMPs 
arc formed, and mpi is the Planck mass. The ratio of 
number abundance of SIMP to that of baryons is then 
Nx/ub - 0.5 X 10-^, where ni, is the number density of 
baryons. 

Diaz-Cruz et al. [iBl have adopted models in which 
the lighter stop ti is the NLSP and the gravitino is the 
LSP. They studied cosmological behaviors of the stop, 
and find that stops are finally hadronized into mesinos 
or sbaryons. They used data on energy levels of normal 
hadron and the estimation of annihilation rate by Kang 
et al. [l4| . and showed an abundance evolution of exotic 
hadrons. 

Jacoby & Nussinov [l7 
relic abundance in Ref. 



have refined the estimate of 
U They assumed that MCPs 
are quarklike, and considered a reaction between SIMPs 
and normal hadrons, i.e., mesons and nucleons, which 
transform constituents or energy levels of exotic hadrons 
X. They suggested that the dissociation of bound states 
of two MCPs must be considered. They roughly took all 
those effects into account, and concluded that the relic 
abundance after the freeze-out of the annihilation at T 
Tc is as small as suggested in Ref. |14 1 . 

Recently the "Quirk Model" has been suggested by 
M. Luty [131 ■ In the model, colored particles transform 
under an additional non-Abelian gauge group. The phe- 
nomenology relating to a new exotic quanta of the cor- 
responding gauge field originally named theton had been 
studied about three decades ago [l^ [13]. They have 
discussed results for the situation in which the mass of 
the lightest quark was heavier than the QCD parame- 
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Taking account of the recent 
an estimate of relic abundances has 
derived, and collider phenomenology of such particles or 
the possibility of detection at LHC or Tevatron has been 
studied [13, [Hj . Strassler and Zurek (2^ have also stud- 
ied the phenomenology of a new confining gauge group 
added to SM kind of hidden- valley models. 

In this paper, we construct a toy model which roughly 
includes effects of all possible processes operating in the 
annihilation at T < Tc. Abundances of SIMPs during 
the annihilation epoch are calculated for four cases of 
different properties of resonant or barely-bound states of 
two-MCP systems and binding energies of SIMPs. We 
show results for each case with different initial abun- 
dances and masses of MCPs. We then suggest that fi- 
nal abundances can be much higher but cannot be much 
smaller than the naive estimations of Refs. [3, [l3l- 

In Sec. ini we describe input physics and assumptions 
adopted in this paper. We show a hypothetical quantum 
mechanical picture of annihilation of SIMPs (Sec. HIK). 
a toy model of the effective cross section for the annihi- 



lation (Sec. ini3), a comparison of the formation rate of 
bound states involving two MCPs with the cosmic ex- 
pansion rate (Sec. |TT]C), the transformation of mesonlike 
SIMPs to nucleonlike SIMPs (Sec. [lip), and the annihi- 
lation cross section at low temperature after the trans- 
formation (Sec. HIE). In Sec. IIIIl we describe models of 
different decay properties of resonances of two MCPs, 
and different binding energies of SIMPs. Effective an- 
nihilation cross sections for the models are derived. In 
Sec. IIV[ we show result of calculations of SIMP abun- 
dance and effects of decay properties of resonances and 
binding energies. In Sec. |Vl we briefly mention the case 
in which heavy hadrons including two MCPs exist sta- 
bly. In Sec. IVIl this work is summarized, and a possible 
impact of results about relic abundance of SIMPs on big 
bang nucleosynthesis is discussed. 

In this paper, the following units are adopted: the 
Planck's constant, the Boltzmann's constant, and the 
light speed are unity, i.e., ?l = 1, fcs = 1 and c = 1. 



II. INPUT PHYSICS 

A. A quantum mechanical picture of resonant (YY) 
production 

Massive colored particles (MCPs), i.e., Y, would be 
confined in colorless states, i.e.. A, at Tc ~ 200 MeV p^ . 
Such states can react with each other and transform into 
a bound state via reactions of the type 



Xi+X2^ (YY) + {7 or qq or qqq ...}, 



(2) 



where Xi and X2 arc states including one Y or Y par- 
ticle, (YY) is a bound state including one Y and one Y 
particle. 

The Schrodinger equation for a two-body wave func- 
tion, xiir) = ripi{r) is given by 



2n dr^ 



1(1 + I) 
2/ir^ 



+ V{r) - E 



Xi{r)=0, (3) 



where fi is the reduced mass, r is the radius of the system, 
V{r) is the central potential at r, E is the kinetic energy, 
and '0(^) is the wave function at r. The centrifugal po- 
tential for an initial state {X1+X2) is small because of its 
large reduced mass. The penetration factor can, there- 
fore, be large even for high angular momentum I. 

Kang et al. [ij] suggested that cross sections for 
Al -I- X2 reaction [Eq. ^] would be as large as that for 
reactions of normal hadrons with their typical values of 
^ [7r(l GeV^^)^]. We note that the cross sections can be 
smaller than their assumption according to a well-known 
concept from nuclear reactions as follows: 

Low- lying resonances or excited states of {QQ)* with 
heavy quarks Q, such as cc or hh, marginally bound with 
respect to the Qq+Qq channel have large width generally. 
The reaction for exotic MCPs is then also considered to 
proceed through such states, i.e., (YY)* . 
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The reaction [Eq. ([2])] for normal nucleons, i.e., A'', is 
Ni + N2 — > (NN) + TT. This is an important reaction 
of pion production and its cross section is measured [l^l . 
The cross section is theoreticahy explained assuming that 
the reaction has a contribution from A(1232) resonance 
(N*) composed of a pion and a nucleon and the existence 
of state NN* HE], [26|. Such a situation is assumed in the 
reaction [Eq. (2)] although the reaction mechanism is 
expected to be much more complicated by many partial 
waves playing roles and many resonant or bound energy 
levels of (VY) as described below. 

A nonresonant component of the reaction might signif- 
icantly contribute to the total cross section. In this study, 
however, we assume that the reaction operates domi- 
nantly through resonances [i^. The cross section would 
then be described using the Breit-Wigner formula (27| . 



where a is the fine structure constant, and AE^ is a typ- 
ical interval between energy levels of resonances. In the 
above equation, A^off ^ is an effective number of 

final states to which the resonance i can transit by emis- 
sion of electric dipole photon. 

A decay width for an emission of Xi is contributed 
possibly from multiple partial waves with large penetra- 



tion probabilities, i.e., Pi.j^i ^ 1 for 
suppressed centrifugal potentials. 
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where E is the kinetic energy in the entrance channel. 
The sum is taken over resonances i. uji = gi/ (3x19x2) is 
the statistical weight factor with ga — 2 -I- 1 the spin 
(Jq) degeneracy factors for a ~ i, i.e., the resonant or 
excited state, heavy hadrons Xi and X2, respectively. 
The total decay width of the state i is given by 



= r 



i,L 



(5) 



where Tij is the partial width for decay of i to j. It is 
given by Ti j = with Ti j i the partial width for 

decay of i to j of relative angular momentum I. Ti^L ^ 



- 2,7 



i,L 

+ Ti^TTTT + ^i,Tv + ■■■ is the decay width of i into all 
light particles (L) composed of SM particles. L does not 
includes Xi and X2. 

In nuclear reactions, the partial width for emissions of 
nucleon or nuclei is expressed by 
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where is the velocity of the exit channel j in reaction 
through the state i, Rij is hadronic interaction radius, 
Pi^j.i is the penetration factor for particles {j) of relative 
angular momentum I. The quantity Of j i is the dimen- 
sionless reduced width which means the probability that 
the state i is described by the state of j of angular mo- 
mentum I. The reduced width is empirically known to 
take values of 0.01 < Oi^j^i < 1 for nuclear states i [27[ . 

There are lots of decay modes of hadrons in which light 
particles are emitted. As we describe hercinbclow, how- 
ever, existing experimental data on heavy quarkonia [28j 
show that the radiative decay width is predominant over 
others. We, therefore, assume that the width for decay 
into all light particles is given by only that of radiative 
decay, i.e., Fj^i = F^^^ in this paper. The width for a 
dipole photon emission is roughly given by 



(7) 



Wc define the sum of reduced widths as 6*1 = J2i Of x i- 
It should be noted that the conservation of angular mo- 
mentum limits possible angular momenta I between Xi 
and X2 to be reached from a resonance of a certain spin 
number. We then obtain an equation, i.e.. 



^^'^ Ri,x\ my 
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(10) 



Experiments of heavy quarkonia [281 ] show that a typ- 
ical decay width of such states are ^ 0(10 MeV). It 
is assumed that intervals between energy levels of res- 
onances lying above the Xi + X2 threshold [or bound 
states of (YY)* lying slightly below the threshold] is 
AEr - 50 - 100 MeV [Tj. These intervals are not so 
large compared with the decay width and the related 
energy scale, i.e., Tq. We assume the equation, i.e., 
\E — Ei\ < Ai?r ^ Fi. The following approximate equa- 
tion then holds: 
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where Ny-^s is the effective number of resonances through 
which the reaction proceeds. In the last equality, we de- 
fined |t 



W7W 



I as an average quantity per resonance 

without distinguishing magnetic substates. 

We assume that spins of Xi and X2 are zero for sim- 
plicity expecting that they are small. A number of reso- 
nances with relative angular momentum I is equal to that 
of magnetic substates, i.e., w ~ (2/ -I- 1) approximately. 
In addition, we assume that typically one resonance for 
each value of I works as an intermediate state in the re- 
action [Eq. ([2])]. We then rewrite the expression for the 
cross section as given by 
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We defined 



as an average quantity per magnetic 



substate. The factor, i.e 
fined as the value ai E ~ 
We used an equation, i.e. 



WJW 



was de- 



E=Ec 

Ec (corresponding to T = Tc). 



Fx Fx 



(r/2) 



- U^yS/Ec, where 

the scahng derives from the width for an emission 
of X [Eq. (Uni)]. 

This equation indicates that the cross section can be 
as large as that of normal hadrons if the following two 
conditions are satisfied: 1. a large number of partial 
waves (Zmax ~ 30 can contribute to the reaction, 

and 2. decay widths to channels for emission of light 
particles L are as large as those for Xi + X2 channels. 

The factor, i.e., /w can be approximately unity at max- 
imum. It can, however, be naturally small if Tx > or 
Fx < r^. The ratio between the decay widths is given 
[Eqs. © and ^] by 
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(13) 



In deriving the first line of this equation, we assumed 
that both of interaction radii at decay of i to X and to 
7 are roughly equal to i?h ^ 1 GeV~^, i.e., the radius of 
heavy hadrons, i.e., Xi. Unless a condition Fx ~ T^ is 
satisfied, values of the cross section reduces by the factor 
/w [see Eq. dm)]. 



B. Effective cross section of YY annihilation 

At temperature T ^ 200 MeV, the annihilation pro- 
ceeds in three stages: a formation of bound states (YY), 
their transitions to lower energy levels, and the annihi- 
lation of YY pair inside the bound states. If a bound 
state tends to be destroyed into separated heavy parti- 
cles at collisions with thermal particles before processed 
by YY annihilation, the state is not an effective path 
for the annihilation [3, [l3l ■ We then introduce an effec- 
tive cross section of annihilation contributed from only 
effective paths for annihilation. It is given by 



(my,T) 



2iiEl!'^Ey^ 
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4.1/w X 10-4 fm^ [Cjmy,r) + l] ' 

-1/2 , , _i 



E, 
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300 MeV 
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200 MeV 



niY 
1 TeV 



-1/2 



(14) 



where l'^^^{mY ,T) is the maximum of angular momenta 
of partial waves available in effectively annihilating YY . 
The maximum angular momentum is given in Sec. IIIIl 
The rate for annihilation of Y is given by 



Tfo 



- dcs 



L for 
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Tea 



(if Teas ^ Tdos) 

(if Fdos > ^ cas } 



(15) 



where Ffor is the rate for bound state (YY) formation, 
and Teas and Fdes are the rates of the bound states for 
cascade down energy levels and for destruction by ther- 
mal background particles, respectively. 

We make a toy model: An annihilation cross section 
for a partial wave I is given by a cross section of bound 
state formation if Teas > Tdes, and is zero if Teas < Tdos- 
The total cross section is then given by Eq. (fT4| with 
ZJj^^(my , T) changing, like the step function, at tempera- 
tures where Teas = Tdes is satisfied for respective angular 
momenta. 



1. Cascade and destruction rates (low-lying I — states) 

The annihilation rate for an ^ = bound state is 
much larger than the destruction rate, i.e.. Teas 2> Fdos- 
Tightly-bound states with small radii have energy eigen- 
values and wave functions determined predominantly by 
the attractive Coulomb-type QCD potential at small 
radii. Their structures are, therefore, very similar to that 
of atoms. 

The annihilation rate for {YY)i^q qq,--- is es- 
timated [iJI similarly to that of positronium. It is 
Tann ^cdM = 5 X 10-^ GeV(aQCD/0.1)5(m/l TeV) 
with aqcD the strong coupling constant. 

The rate for destruction by thermal photons (YY)i^q + 
7 — >■ Xi -|- X2 is estimated using the detailed balance re- 
lation [2^ as follows: The cross section for the inverse re- 
action, i.e., X1+X2 (F?);=o+7, through an interme- 
diate state (YY)* with angular momentum I is roughly 
given by 



^(2/ + l)/w 
myE^^^E^/^ ' 



(16) 



Supposing the usual atomic model, the number of 
states (YY)* of angular momentum I is ^ /,„ax — ^ + 1, 
where /,nax is the maximum angular momentum of all 
bound states. There are thus many exit channels gen- 
erally. Assuming that all partial cross sections for final 
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states of a given I are equal, the cross section for a fixed 
final state without distinguishing substates is given by 



TT {21 + 1) /« 



(17) 



The destruction rate is derived using this cross section 
and the detailed balance relation. It is then given by 
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1 7-3/2 
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3 X 10^" GeV 
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exp{-Eth/T) 
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300 MeV 
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200 MeV 
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4 X 10- 



(18) 



where is the background photon number density, {ac) 
is the thermal average value of ac, a{E) is the cross sec- 
tion at energy _E, i?th = Q is the energy threshold of the 
destruction reaction, and Q > is the reaction Q-value. 

The binding energy of states with main quan- 
tum number n is E^ ~ /iaQCD/(2n^) = 
2.5 GeV/n2(aQCD/0.1)2(TOy/l TeV). It can be 
significantly higher than the confinement temperature 
scale. In that case, the destruction rate is hindered 
since only photons energetic enough to dissociate the 
state can react. The bound I = states produced via 
Xi -|- X2 collisions thus annihilate before being destroyed 
by photons. 



2. Cascade and destruction rates (high-lymg I > 1 states) 

The annihilation of Z > 1 states is hindered because of 
large angular momenta and low values of wave functions 
at r = [4^ . The fate of the states then depends on the 
rate for transition to states of lower energy and angular 
momentum, and that for destruction by photons. 

Because of color confinement, the predominant bound- 
bound transition by the strong interaction associates 
with transitions of at least second order in perturbative 
color electric or magnetic 2'-pole moments (El or Ml, re- 
spectively) [lOl- The final state is a bound state, i.e., 
(YY)* , and light hadrons. Transitions of A/ > 1 to final 
states with smaller angular momenta are hindered since 
they need high multipole transitions. 

An effective process reducing angular momenta is 
the dipole photon emission if the Y particle has a fi- 
nite charge. The transition rate is given by Teas ~ 
arg^j^j.Ai?'^ ^ TOyaaQQp/lG. An only difference from 
atomic dipole transition rate originates from the tighter 
binding by the Coulomb-like QCD force. The QCD 
force shrinks wave functions and enhances the transition 
strength. 



In particle data [2^, it is seen that bottomonia (char- 
monia) bound against decay into BB {DD) channels pre- 
dominantly decay into bottomonia (charmonia) of lower 
levels by emitting one photon or two pions except for a 
decay of 6 or 6 (c or c) itself. The relatively strong tran- 
sition associated with a two-pion emission between two 
color-singlet states is via color El-El interaction and the 
combined angular momentum of pions is ^ = [3fl] . The 
most strong transition with a change in angular momen- 
tum is, therefore, via a photon emission of mainly 1 = 1. 
We then assume that the most effective process reduc- 
ing angular momenta of the system is the dipole photon 
emission, and consider this process only. 

In addition to the spontaneous dipole photon emis- 
sion, collisional deexcitation reactions can be important 
[l3,[3- Jacoby & Nussinov estimated that the cross sec- 
tion for collisional deexcitation is smaller than that for 
breakup of {YY) into Xi -f X2 by - lO'^-lO-^. The 
deexcitation reaction triggered by a photon, i.e.. 



{YY) + j^{YYy + 



(19) 



through resonant states is expected to be hindered by a 
factor /; = [r^r;/(r/2)2] [see Eq. dH])]. For example, 
when the radiative decay width is typically 10^^ times 
the total width, the hindrance factor is /4 ~ 10~^. The 
cross section is then lO^'^ times as large as that for the 
{YY) breakup [the inverse reaction of Eq. 

The resonant deexcitation reaction triggered by a pion, 
i.e.. 



{YY)+TT^ {YY)' + TT, 



(20) 



is hindered. This is because the pion in the final state 
has I = 2 relative angular momentum (l7l |. Reactions 
involving small reduced masses and finite angular mo- 
menta are, however, suppressed by a large Centrifugal 
potential [cf. Eq. ([3])]. Thus, the resonant deexcitation 
cross sections would be smaller than that for the {YY) 
breakup. 

In this study we assume that radiative resonant re- 
actions have a strength of = lO^^/w, and are the 
dominant deexcitation process as one example case. The 
deexcitation cross section is, however, uncertain signifi- 
cantly. We then adopt the rate given [cf. Eq. ([T8| ] by 
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(21) 



As for destruction processes of / > 1 bound states, we 
assume that bound-bound transitions are mainly through 
the dipole photon absorption which is inverse of the emis- 
sion process mentioned above. Bound-free transitions 
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would be through absorptions of photons with Z = 1 or a 
two-pion pair with 1 = 0. The latter is neglected in the 
present estimation. The transition rate is then given by 
Eq. (HH). 



C. Comparison of {YY) formation rate and the 
Hubble expansion rate 



The YY annihilation does not operate effectively if the 
formation rate is smaller than the expansion rate, i.e., 
H > Ffor- In general, bound states with smaller binding 
energies become less subject to destruction by thermal 
particles at lower temperatures. In such low temperature 
environments, however, the bound state formation might 
be already inefficient because the ratio of rates scales as 



In this subsection, we describe a way to estimate the 
relic abundances of SIMP, i.e., X, for several possible sit- 
uations. In what follows, the MCP, i.e., Y, is assumed to 
be quarklike as an example. The most important process 
of bound state formation and its Q- value are then 

Yq + Yq^{YY) + {j, qq, ...}, 
Q = {Eb[{YY)]+Eb{j, qq, ...)} - [EB{Yq) + EsiYq)] 
^ Eb[{YY)] - or or + 2Am(rg), (22) 

where E-sii) is the binding energy of species i. The bind- 
ing energy is the energy level of (compound) particle mea- 
sured from the sum of rest masses of constituents. The 
gain in mass by a binding can be defined as Am(i) = 
—Esii), and we have assumed Am{Yq) = Am{Yq). 
From Eq. (|14p , the reaction rate is given by 



v) 



(23) 



We note that, if Q-values were negative for production 
of some (YY)* states, then their rates are small since 
only photons or pions energetic enough to induce the 
endoergic reaction can play a role. 

Under the assumption of /w ^ 0.01 and Z^^^ ~ 30, the 
reaction rate of bound state formation is given by 
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112C(3) 3 ny 

A _9 
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= 8 X 10"^^ GeV 



300 MeV 
T 

40 MeV 



0.01 

1/2 



my 
1 TeV 

77 



-3/2 



^off 



6 X 10~io 
2 



10-6 



1 



30 



(24) 



where (11/4) is a factor concerning annihilation at 
T < me/3 with me the mass of electron, ^(3) = 1.2021 
is the third Riemann's zeta function, 77 is the baryon to 
photon number ratio of the present Universe, ny and rib 
are the number densities of Y and baryon, respectively. 
The Hubble expansion rate is given by 



H = 



2 1/2 2 



3/5 m-p\ 



7.1 X 10" 



GeV 



V10.75/ V 



T 



40 MeV 



(25) 



D. Transition of Yq Yqq 

At T '-^ 40 GeV, the coUisional annihilation of thermal 
baryons containing quarks with antibaryons is thought 
to freeze out. Its number density is then fixed except for 
the dilution by cosmic expansion. Abundances of antin- 
ucleons keep decreasing via annihilation with abundant 
nucleons. This is because there are extra quarks over an- 
tiquarks, that is a finite baryon number density in the 
Universe. The reaction 

Yq + qqq O Yqq + qq (26) 

then drives the Yq hadron to Yqq [H [H, 113 ■ The tem- 
perature for operation of this transition is estimated us- 
ing the Saha equation, i.e.. 



fly qTlqqq 



gVqgqqq ( fJ-Yq+qqq 



nyqqUqq QYqqQqq \f^Yqq+qq 



3/2 



exp(-Q/T), (27) 



where Ui is the number density of species i. 

Only pions and nucleons are considered as mesons and 
baryons, respectively, in a thermal bath. The equations, 
qq = TT, Jqq = 0, qqq = N and Jqqq = 1/2 then hold. 
Heavy hadrons, i.e., Yqq and Yq, are assumed to exist 
exclusively in only the ground state. Using the rough 
assumption that the number abundance of pion is given 
by the classical Boltzmann distribution, we obtain a ratio 
between number abundances of Yqq and Yq: 



nyqq 
nyq 



nN 1 / ™^ 

7l7r 2 WTljv 



3/2 



cxp(Q/T) 



3v/2^ 



/ rp \3/2 


Q + mT, 


exp 


T 


\mN J 



V 



6.2 X 10- 



T 



37 MeV 



3/2 



f exp [(Q + m^)/T] \ 
^ \exp(mAr/37 MeV) J 



(28) 



where mpf = 0.94 GeV is the nucleon mass, and m^r = 
0.14 GeV is the pion mass. 

In this study, we take two different Q- values, i.e., 
Q = [EniYqq) + Enim)] " [EsiYq) + E^iqqq)]. The 
Saha equation [Eq. p7|) ] describes that Yq particles are 
processed into Yqq particles at Ttra ^ 40 MeV depending 
on the Q-value. Ttra is the transition temperature below 
which the condition, i.e., nyqq/nyq > 1, is satisfied. An 
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important point is that the Q-value is large compared 
with the temperature T . If a Q- value were smaller, the 
Y particle is mainly confined in Yq until the temperature 
becomes much smaller (see footnote 1 on page 2 of Ref. 
0). 

Since there are only small abundances of antibaryons 
at this low temperature, the Y particles are still con- 
fined in Yq. They do not experience the corresponding 
reaction triggered by antibaryons. 



E. Annihilation through Yqq + Yq reaction 

After the transformation, i.e., Yq — > Yqq, the main 
reaction for bound state formation changes to 



Yqq + Yq'^ [YY] + qqq. 



(29) 



This forward reaction, however, proceeds with a nucleon 
creation. Therefore, unless energy levels of Yqq and Yq 
are rather high or the binding energy of (YY) is enough 
large, the inverse reaction tends to dominate because of 
negative Q- value. 

The reaction rate of bound state formation is estimated 
as follows similarly to the calculation of Eq. (p4)) . As 
the temperature decreases and the number density of 
energetic X particles decreases, resonant reaction rates 
for resonances of high energies are expected to become 
smaller. Such an inactivation of resonant process is, how- 
ever, very complicated and not studied in the present 
study. 



1. Case of [Am{Yqq) + Am(yg)] > mjv 



The rate is given by 



GOV (^)- 

Z7 \ -1/2 / 



300 MeV 



30 



6 X 10-1" ) Vl-7 X lO-'^ 

(30) 



20 McV / ' 



where we roughly assumed that reactions from all initial 
partial waves contribute to the {YY) formation. Destruc- 
tions of {YY) states after their formation by nucleons or 
photons were, therefore, neglected since the temperature 
is already rather low. 

Final abundances of Yqq and Yq are estimated as 
values which correspond to the same reaction rates as 
the Hubble expansion rate. The ratio between the Y 
abundance and that of baryon is thus estimated to be 
ny/nij ^ 2/~^ X 10^^. The bound state formation after 
the transformation of Yq — > Yqq is thus somewhat inef- 
ficient in reducing the Y abundance. Relic abundances 



of exotic hadrons including one Y or Y particle then de- 
crease affected by only cosmic expansion since the time 
of transformation to Yqq. 



2. Case of [l\m{Yqq) + Am{Yq)\ < tun 

The reaction rate is estimated to be smaller than in the 
previous case because of smaller (or negative) Q-values. 

We define e';^^ and a'^ as the threshold energy and the 
partial cross section, respectively, of reaction [Eq. (|29l) ] 
whose final state is a (YY) state of angular momentum 
I and a qqq (nucleon). The reaction rate is then given 
by the sum of all angular momenta of (YY) in the final 
state, i.e.. 



Tfor 



ny- 



TO7^E(2^ + 1) /I V^exp(-x)d. 



ny 



Ef. I 

3/2 j-,1/2 

I 



ij- E' ,_, . Je^^ T 



7/w X 10-22 Gcv /: 



2 

mill 



my \-3/2 



VI TeV/ 



Ec 



300 MeV 



-1/2 



6 X 10-10 



^ ' 7 X 10-5 ) Uo McV / '■ 



(31) 



where lm\i\ is the minimum angular momentum of {YY) 
which have negative Q- values, so that reactions produc- 
ing such states are endoergic. First and second term in 
square bracket in the third line correspond to the cross 
sections for partial waves with I < /min and I > Imin, 
respectively. The former component has zero threshold 
energies so that it does not suffer from Boltzmann sup- 
pression factors relating to a high threshold energy. The 
second term is neglected in the fourth line since it is gen- 
erally subdominant to the first term, [s^ 

Final abundances of Yqq and Yq are estimated to be 
ny/nt ^ 7/-^?-?^ X 10"^. This value is significantly 
larger than that in the previous case of [Am{Yqq) + 
Am{Yq)] > mpf. 



III. MODEL 

We describe a model for an estimation of the relic 
abundance of the exotic heavy hadron X. Table U lists 
parameters for four cases we consider in this study. Case 
1 is our standard case. We chose them in order to clearly 
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TABLE I: Parameters for four cases. 





/w 


Am(Fg) 


Am{Yqq) 


Case 1 


0.01 


0.3 GeV 


0.7 GeV 


Case 2 


1 


0.3 GeV 


0.7 GeV 


Case 3 


0.01 








Case 4 


1 









SIMPs (i?-hadrons) in colliders were studied. For coni- 
parison, in Cases 3 and 4 the mass gains are set to be 
zero supposing that their amplitudes are not large. 

We note that the Q-values of the reaction [Eq. (|26p ] 
are Q = rriN — — 0.4 GeV in Cases 1 and 2, and Q = 
rriN'-mT^ in Cases 3 and 4, respectively. Such a difference 
in Q-valucs leads to different transition temperatures [cf. 
Eq. (UHl)], i.e., Ttra 20 MeV (Cases 1 and 2), and Tt^ ~ 
37 McV (Cases 3 and 4). 



show effects of resonance cross section and energy levels 
of the X. 



C. Effective cross section for Yq + Yq 



A. Resonance factor /„ 

Values for resonance factor arc adopted as follows. 
The particle data on charmonia and bottomonia psj pro- 
vide a possibly useful information. We consider reactions 
relating to heavy c or 6 quarks of the type of Eq. 
The lowest resonance of charmonium above the levels of 
two-ZJ-mesons states {D^ = cd, D'^ = cu) is V'(3770). Its 
fraction of partial decay width for D + Z) channel relative 
to the total width is (93^g) x 10~^, while the sum of those 
for other cc mesons + light particle channels (7, tttt, r/, 7r°) 
is 1.5 X lO^'^ 1^. The factor of resonance decay widths 

is then /„ ~ r^(3770),r>r^(3770),{7,7r7r,...}/(r.v>(3770)/2) = 

0.06. The rate for the reaction of type [Eq. ([2])] through 
this resonance is, therefore, smaller by the factor /w than 
the naive estimate of the cross section corresponding to 
the hadronic size. 

The lowest resonance of bottomonium above the levels 
of two-_B-mesons states (i?+ = ub, = db) is T(4S') [or 
T(10580)]. Its fraction of partial decay width for B + B 
channel relative to the total width is larger than 0.96, 
while the sum of those for other bb mesons + light particle 
channels (tttt, rj) is smaller than 4.4 x 10""^ [2^. The fac- 
tor is then = rx(4S).SrT(4S),{7r7r,,,...}/(rT(4S)/2)^ < 

0.018. 

As for reactions of two bound states involving Y par- 
ticles, resonance widths have large uncertainties. We as- 
sume that the factor is /w — 0.01 taking account of the 
value, i.e., /w < 0.018, derived from particle data on the 
bb resonant state. This state is composed of the most 
heavy quark b whose data on resonance decay widths are 
available. The factor for X hadrons may be larger or 
smaller than /w = 0.01 in fact. We, therefore, also show 
results for the case of = 1 as well as those of = 0.01. 



B. Energy levels of Yq and Yqq 

We take two different reaction Q-values in order to 
show that Q- values or energy levels of exotic hadrons, i.e., 
Yq and Yqq and so on, are significant in determining final 
abundance of SIMPs. In Cases 1 and 2, we take the mass 
gains of Am{Yq) = 0.3 GeV and Am{Yqq) = 0.7 GeV as 
assumed in Ref. [3l| in which searches for scatterings of 



We treat the annihilation through the Yq + Yq re- 
action in the temperature range of T > Ttra- The ef- 
fective cross sections [Eq. p4)) ] for Cases 1-4 are given 
using l'^^^^{mY ,T) determined by a procedure described 
in Sec. HlB. Here, each l'^^^^{mY,T) is determined with 
values of critical temperatures, i.e., Tcri(l). These critical 
temperatures are defined by points at which the condi- 
tion, i.e., Fcas = Pdcs, is satisfied. Whether the cascade 
or destruction dominates then changes at Tc,i(0 for each 
angular momentum I. A cross section for partial wave I 
is maximal if Fcas > Fdos, and is zero if Fcas < Fdes- 

In addition, we adopt an ansatz on the quantity 
ZJ^^(my,T) at temperatures below the critical value 
7cri(7) for Cases 1, 2 and 3. The ansatz is a linear func- 
tion of T, i.e.. 



=7- 

max ' 



rcri(7) - T 



(32) 



where Zmax = 30 is fixed in this study. We calculate 
critical temperatures Tcri(0 below. 



1. Comparison of rates for cascade and destruction 

The destruction rate is given by Eq. ([T5)) . The thresh- 
old energy is Eth ^ 2.5 GcY / n"^ {my / 1 TeV) -I- 0.6 GeV 
(Cases 1 and 2) and Eth = 2.5 GeV/n2(my/l TeV) 
(Cases 3 and 4). 

The cascade rate is given by sums of spontaneous rates, 
i.e, F^Pg, [Eqs. (f33| and (|34l) ] and colhsional rates, i.e., 

r™ , [Eq. m]- 

For the spontaneous cascade rate of states with small 
I, we roughly take the rate for a transition from state 
(n, / = n — 1) to {n — l,l — 1) as a typical rate for dipole 
emission from a state with angular momentum I. It is 
described by 



psp 



Q![rBohr('^ 



!)]■' [Ejiin) - Ejiin - 1) 
(Z + 1/2)3 



(Z-hl)6P 



(33) 



where rBohrin) and £'b(?^) are the Bohr radius and the 
energy eigenvalue of state with main quantum number n. 

States of I ~ ^max with hadronic sizes in space and 
small binding energies, on the other hand, take longer 
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i . t < . . \ 

50 100 150 200 250 300 

T (MeV) 



FIG. 1: Rates for cascading down energy levels, i.e., Teas 
(solid lines), and for destruction by thermal photon, i.e., Fdca 
(dashed), of bound states with / ~ n — 1, as a function of 
temperature. The factors of resonance widths are chosen to 
be /„ = 0.01 and /4 = 10"''. The mass gain of Y q hadron 
is Am(yq) = 0.3 GeV. Masses of the Y particle are taken 
to be 500 GeV (thin lines), 1 TeV (intermediate) and 3 TeV 
(thick). Angular momenta of {YY) are taken to be 1, 5 and 
^max (from the top line to the bottom for Teas, and the bottom 
to the top for Fdos). 



time to cascade down [TJ: 113 • The rate for an electrically 
charged Y case is roughly given by 

psp ^ "^^had 
cas ^ 1/2 2 

^ 10-«Gevf^Vf^)"'^'f^)"'. 

VI GcV y V 0.1 / VI TcV/ ' 

(34) 

where Ahad is the energy scale of hadronic interaction. 

Figure [1] shows Teas (sohd lines) and Fdcs (dashed) of 
bound states with I ^ n — \, 0jS &- function of tempera- 
ture. The colored Y particle has been assumed to have 
an electric charge. This figure is for Case 1: The fac- 
tors of resonance widths are /w = 0.01 and = 10""*, 
and the mass gain of Yq hadron is Am(yg) = 0.3 GeV. 
Masses of the Y particle are taken to be 500 GeV (thin 
lines), 1 TeV (intermediate) and 3 TeV (thick). Angular 
momenta of {YY^ are taken to be 1, 5 and ^max (f^rom 
the top line to the bottom for Teas, and the bottom to 
the top for Fdcs). 

The I = 1 excited states produced by the hadronic 
reaction [Eq. ([2|)] arc found to instantaneously trans- 
form to / = ground state in a temperature range of 
T < 300 MeV, as seen from Fcas > Fdcs- The / = 5 states 
also tend to transform to Z = 4 state rather than dissoci- 
ated in a photo-reaction at T < 141 — 186. Critical tem- 
peratures change slightly depending upon the mass of the 
colored particle. On the other hand, as for (YY) states 
barely bound relative to the Yq + Yq separation channel. 




T (MeV) 



FIG. 2: Same as in Fig. [T]for the case of Am{Yq) = 0. 



they can have large angular momenta of I ^ Imax and 
small positive reaction Q-values of Q ^ 100 MeV. Be- 
cause of the fragility of such barely bound states against 
background photons, they tend to be destroyed before 
cascading down to lower energy states until the temper- 
ature decreases to T ~ 0(10) MeV. 

Figure [5] shows Fcas (solid lines) and Fdcs (dashed) sim- 
ilar to Fig. [l]for Case 3: The factors of resonance widths 
are /w = 0.01 and /4 = 10"^, and the mass gain of Yq 
hadron is Am{Yq) = 0. The smaller value of ATn{Yq) 
leads to smaller Q-values for respective bound states. 
The bound states are, therefore, more easily destroyed 
by background photons. Destruction rates are thus gen- 
erally larger than in Case 1. The I = 1 excited states 
instantaneously transform to the I = ground state sim- 
ilarly to Case 1. The I ~ 5 states tend to experience 
transitions to lower states at T < 11 — 54. Barely bound 
(YY) states tend to be dissociated by background pho- 
tons until the temperature decreases to T ^ 0(10) MeV. 



2. Estimation of effective cross sections 

Table HIl shows critical temperatures in Case 1. The 
first column shows the mass of the Y particle. The second 
to ninth columns list critical temperatures for respective 
I values. Energy levels naively based on the Coulomb 
models were assumed. States with main quantum num- 
bers n > 8 would be mainly affected by not a Coulomb- 
type QCD potential but a linear potential. We do not 
treat details on loosely-bound energy states originating 
from this linear potential. The energy levels and wave 
functions, however, should be studied for various likely 
potential models in future to calculate formation cross 
sections and rates of transitions to lower energy states 
precisely. 

Table Hill shows critical temperatures in Case 2 similar 
to Table HIl In Case 2, the setup is the same as in Case 
1 except for the difference in the resonance factor. The 



10 



TABLE II: Critical temperature of annihilation for Case 1. 

niY (TeV) Tcri (MeV) 

; = 01 2 3 4 5 6 7 

0.5 608 285 199 165 149 141 137 135 

1 887 375 243 190 164 151 144 140 

3 1820 670 384 271 216 186 169 158 



(YY) bound states is small since dissoeiations triggered 
by background photons are weak as in the inverse re- 
action. Small rates of the (YY) bound state formation 
in the case of neutral Y particles would then result in 
more early freeze-out of annihilation in the Universe. Fi- 
nal abundances of exotic heavy hadrons arc rcsultingly 
larger than in the charged case. 



TABLE III: Critical temperature of annihilation for Case 2. 



ruY (TeV) 


Tcri (MeV) 


;=01 2 3 4 5 6 7 


0.5 
1 
3 


231 165 148 142 138 136 135 134 
318 199 166 153 146 141 139 137 
646 331 237 197 176 164 156 150 


TABLE IV: Critical temperature of annihilation for Case 3. 


ruY (TeV) 


Tcri (MeV) 

/=01 2 3 4567 


0.5 
1 
3 


381 110 49 27 17 11 8 6 
659 195 88 48 30 20 15 11 
1599 490 226 126 79 54 39 29 


TABLE V: Critical temperature of annihilation for Case 4. 


ruY (TeV) 


Tcri (MeV) 

1^0 1 2 3 4 5 6 7 


0.5 
1 
3 


99 37 19 11 8 6 4 3 
187 71 38 23 16 12 9 7 
516 203 109 68 47 34 26 20 



factor in Case 2, i.e., /w = 1, is 100 times larger than that 
in Case 1. Destruction rates Fdes in Case 2 are then 100 
times as large as values in Fig. [TJ Because of relatively 
larger destruction rates, critical temperatures change to 
lower values. 

Table IIVI shows critical temperatures in Case 3. Be- 
cause of reaction Q- values smaller than in Case 1, the 
destruction process by background photons is stronger. 
This leads to delayed transition, or smaller critical tem- 
peratures than in Table HIl 

TablelVlshows critical temperatures in Case 4. Because 
of reaction Q- values smaller than in Case 2, the destruc- 
tion process is stronger. The transition is then delayed 
and critical temperatures are smaller than in Table IIIII 

If the colored Y particle were neutral, its rate for cas- 
cade down would be small and the transition would be 
via emissions of light mesons. Formation cross sections 
of (YY) bound states [Eq. (|12p ] are also small since 
channels of dipolc photon emissions are not open in res- 
onant reactions. In addition, a dissociation rate of the 



IV. RESULT 

Relic abundances of exotic heavy hadrons are calcu- 
lated using the annihilation rates described in Sec. IIIII 
We denote hadrons which include one Y particle by X, 
and hadrons which include one Y particle hy X. It should 
be noted that Y particles are not always antiparticles of 
Ys. The time evolutioir of the number density of heavy 
hadron X {X) is described by the following integrated 
Boltzmann equation [3^ : 



dYx 
dt 



dYx 
dt 



^-.s{a{mY,T)v)YxYx, 



(35) 



where Yi = n.i/s was defined as a measure of the number 
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of species i = X and X per comoving volume 
is the number density of i, and s is the entropy density 
of the Universe, {airny ,T)v) is the effective reaction 
rate for YY annihilation via formation of bound states 
of MCPs (YY) by two SIMPs. In this equation, the 
inverse reaction of annihilation, i.e., the heavy hadron 
production, has been neglected. Since the temperature 
of the Universe T < 200 MeV is much lower than the 
mass of MCP m ^ 0{1 TeV), such a production is never 
operative. 

Effects of the destruction and decxcitation of bound 
states (YY) by thermal particles, and the transition of 
heavy hadron species, i.e., Yq — Yqq, are taken into ac- 
count in the reaction rate {a{mY ,T)v). We assume that 
X and X have the same number density at the temper- 
ature of the QCD phase transition Tq ~ 200 MeV [Slf . 
Then, the equations for time evolutions of Yx and Yx 
arc the same, and the equation, i.e., Yx = Yx, holds. 
The equation is described by 



dY: 



X 



dt 



= ~s{a{mY,T)v)Yx- 



(36) 



This equation is solved as a function of temperature T 



Yx{T) = Yx{Tc) 



l + Yx{Tc)m^u 



Tc 



X / dT^g,s{T){a{mY 



,(37) 



where ^/gMT) = ^54^7 (l + 3^;^%r^) was de- 
fined. 
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100 10 1 

Temperature T [MeV] 

FIG. 3; Yields of heavy hadrons X and X as a function of 
temperature in the region T = 200 MeV-1 MeV for Case 1. 
The mass of the Y is fixed to be my = 1 TeV. The sohd, 
dashed and dot-dashed hnes correspond to the initial yields 
of heavy hadrons of y(rc) = 10"", 10"" and 10"", respec- 
tively. 
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FIG. 4: Same as in Fig.[3]for Case 2 for my = 1 TeV. 



In Figs. [SHll we plot calculated abundances of X and 
X, i.e., Y{T) = Yx{T) + Y^iT) = 2YxiT) for Cases 1-4 
(see Table m, respectively. We took three initial abun- 
dances of Y{Tc) = 10-^^, 10-1-* and 10"l^ and results 
are shown by the solid, dashed, and dot-dashed lines, re- 
spectively. Fig. 13] shows results for our standard case, 
i.e.. Case 1. 

In Fig. m the value is 100 times lager than that in 
Case 1, and the mass gains of Yq and Yqq are the same 
as in Case 1. This difference results in relic abundances 
smaller than in Fig. [3] by a factor of ~100. In Fig. 21 the 
final abundances of exotic hadrons do not depend at all 
on the initial values of Y{Tq). 

In Fig. [5j while the /w value is the same as in Fig. [3l 
the mass gains of exotic hadrons are assumed to be zero. 
Because of the zero mass gains, the bound state forma- 
tion rates arc smaller and the resulting relic abundances 




100 10 1 

Temperature T [MeV] 

FIG. 5: Same as in Fig. [3] for Case 3 for my — 1 TeV. 
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FIG. 6: Same as in Fig. [3] for Case 4 for my = 1 TeV. 



are larger than those in Case 1. The YY annihilation 
at the hadronization of the MCPs Y and Y is not so 
important to the final abundances in Case 3. 

Finally, Fig. [5] shows results for Case 4 in which the 
/w value is 100 times lager than in Case 1 and the mass 
gains are assumed to be zero. The final abundances in 
Case 4 are similar to those in Case 1. This is because the 
effects of the resonance factor and the mass gains work 
in different directions and they get balanced out. 

Nonsmooth behaviors can be observed in Figs. [5] and 
Fig. [B] They come from discontinuities in adopted 
^max('7iF,T) values which lead to discontinuities in the 
effective cross sections [Eq. (|14|) ] and reaction rates. The 
discontinuity would be an artificial fake made in our sim- 
plified toy model, and not be physical. We expect that 
more smooth cross sections would realize in more dedi- 
cated models. 

In Fig. [51 plateau yields are seen at T ~ 200 MeV 
(leftmost part). This flat structure results from 
{a{mY,T)v) = 0. It reflects that rates for cascading 
down of (YY) states are smaller than those for their de- 
struction by thermal particles for all partial waves at the 



12 



: 1 1 1 1 — 1 — 1 — 1 — 


-1 1 1 1 1 1 — 1 — 1 — 1 1 1 : 

my = O.STeV : 




1 TeV ■ 




3TeV 















100 10 1 

Temperature T [MeV] 



FIG. 7: Yields of heavy hadrons X and X as a function of 
temperature in the region T = 200 MeV-1 MeV for Case 1. 
The initial yield is fixed to be Y{Tc) = 10"". The solid, 
dashed and dot-dashed lines correspond to the masses of Y 
of rriY = 0.5, 1 and 3 TeV, respectively. 



V. EFFECT OF (YVq) HEAVY HADRON 

We note that states involving two exotic heavy colored 
particles, i.e., (YYq), would naturally exist if the YY 
annihilation does not operates inside the hadron. On the 
other hand, if two Y particles interact inside the hadron 
and transform into SM particles, the particles (YYq) play 
a role in reducing the abundance of Y. Productions or 
detections of such hadrons with double Ys in colliders 
might be rather difficult. They are, however, thought to 
have been produced and also destroyed in a hot epoch of 
the early Universe. 

Although we neglect these states in this study, we com- 
ment on their possible effects assuming that they are 
long-lived. The most important thing regarding the ef- 
fects would be reaction Q- values which stem from energy 
levels of relevant hadrons. Possible processes for produc- 
tion of (YYq) [and (YYq)] are the following three (and 
similar three for antihadron): 



TABLE VI: Calculated number ratios of heavy hadrons to 
normal baryons rix+x/ rib at T = 1 MeV for four Cases (see 
Table |l]) and different masses my and initial yields of heavy 
hadrons, i.e., Y(Tc). 



Case 


y(rc) 


my = 0.5 TeV 


1 TeV 




3 TeV 






10-" 


3.4 X IQ-'' 


9.2 X 10" 


6 


4.2 X 10" 


5 


Case 1 


10-" 


3.3 X IQ-'' 


8.6 X 10" 


6 


3.2 X 10" 


5 




10-" 


2.6 X 10-® 


5.2 X 10" 


6 


9.2 X 10" 


6 




10-" 


3.5 X 1Q-* 


9.6 X 10" 


8 


4.5 X 10" 


7 


Case 2 


10-" 


3.5 X 1Q-* 


9.6 X 10" 


8 


4.5 X 10" 


7 




10-" 


3.5 X 10-* 


9.5 X 10" 


8 


4.4 X 10" 


7 




10-" 


3.3 X 10-'' 


3.9 X 10" 


4 


5.7 X 10" 


4 


Case 3 


10-" 


9.3 X IQ-'^ 


9.8 X 10" 


5 


1.1 X 10" 


4 




10-" 


1.1 X 10-5 


1.2 X 10" 


5 


1.2 X 10" 


5 




10-" 


2.2 X 10-5 


2.0 X 10" 


5 


2.1 X 10" 


5 


Case 4 


10-" 


1.9 X 10-5 


1.8 X 10" 


5 


1.8 X 10" 


5 




10-" 


7.7 X 10-" 


7.5 X 10" 


6 


7.6 X 10" 


6 



high temperature. 

For illustration of the mass dependence of the relic 
abundance, in Fig. [7] we plot yields Y{T) for cases of dif- 
ferent masses. The initial yield is fixed to be Y(Tc) = 
-|^g-i4 rpj_^g solid, dashed, and dot-dashed lines corre- 
spond to the masses of Y of my — 0.5, 1 and 3 TcV, 
respectively. One can see that the larger mass results in 
the larger relic abundance. This behavior derives from 
a dependence of the formation rate on mass my- by the 
negative power [cf. Eq. (^5]) ]. 

In Table IIVI we summarize calculated number ratios 
of heavy hadrons to normal baryons, i.e., 
T = 1 MeV. For Cases 1-4, results are shown for different 
masses and initial yields. 



A. Yq + Yq iT>Ttr^) 



Yq + Yq^ (YYq) + qqq, 
Q = EB[{YYq)]-mN + 2Am{Yq). (38) 

In Cases 1 and 2, the Q- value is given by Q ~ Eb (YYq) — 
TOjv + 0.6 GcV. Because of the height of the resonance 
level of Yq, and the mass of nucleon in the exit channel, 
the Q-value is expected to be negative. This reaction is 
then not likely an available path for two Y (or Y) parti- 
cles to be bound. If a Q-value is positive and relatively 
large, e.g. of 0(100 MeV), however, the colored heavy 
particles are effectively processed through this reaction, 
and locked in the (YYq) hadron. 

The (YYq) hadron thus produced can be processed 
through two reactions, i.e., 

(YYq) + Yq o (YYY) + {7, gg, ...}, 
Q - Eb[{YYY)] - m{^, ,5,...} - {EB[{YYq)] - Am{Yq)], 

(39) 

and 

(YYq) + Yq^ (YY) + Yqq, 
Q = {EbUYY)] - Am{Yqq)} - {EniYYq) - Am[{Yq)]}. 

(40) 

Whether these two effectively operate depends on their 
reaction Q-values. However, both reactions possibly have 
positive Q- values. The colored particle Y would then be 
predominantly included into (YYY) which is long-lived 
and (YY) which eventually annihilates. 

The (YYY) hadron is expected to be rather tightly 
bound because of its large mass. The only possible reac- 
tion processing low-lying states of (YYY) with positive 
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Q- value is then 

{YYY) +Yq^ {YY) + {YYq), 
Q = {Eji[{YY)] + E^[{YYq)]} - {Eji[{YYY)] - Am{Yq)}. 

(41) 

Energy levels of the entrance channel including {YYY) 
are typically lower than those of the exit channel includ- 
ing (YY), although somewhat excited {YYY) states may 
be transformed into tightly-bound low-energy levels of 
{YY). Then, a significant fraction of initially existent Y 
might be confined in {YYY) particles. The {YYY) par- 
ticle would escape from reducing its abundance since the 
annihilation {YYY)+{YYY) is hindered by somewhat 
larger mass and significantly smaller size of the targets 
composed of three heavy Ys than those composed of a 
single Y and one light q or two light qs. 

Final abundances of {YYY) and {YYY) in this case 
would be roughly of the order of initial Y abundance 
fixed at freeze-out of perturbative annihilation, i.e., 
'in(YYY)/''^b ^ 10"^'*. Depending on situations, final 
abundances of {YYq) can be high, i.e., 2n(YYq)/nb ^ 
10-". 



B. 



Yqq + Yq (T ~ Ttra) 



Yqq + Yq^ {YYq) + {-l, qq, ...}, 
Q = EniiYYq)] - m^^^ + [Am{Yqq) + Am{Yq)]. 

(42) 

In Cases 1 and 2. the Q- value is given by Q < 
£'b[(^^9)] + 1 GeV depending upon light particles in 
the final state. The Q-value might be positive and large. 
In our estimation, abundances of Yqq and Yq are compa- 
rable around the transition temperature Ttra ~ 20 MeV. 
At this temperature, this reaction can operate. If the 
YY annihilation through {YYq)+Yq or {YYY)+Yq is 
hindered by possible large negative Q-values for the 
reactions, the Y particles would survive being domi- 
nantly locked in the Yq hadron. If the reaction [Eq. 
(|^ ] effectively operates, final abundances of (YYq) or 
(YYY), and Yq would be roughly of the order of Yq 
abundance at the transition temperature or less, i.e., 
3n(^YYY)/nb < lO^MO-"* or 2n(^YYq)/nb < 10-^-10"'*, 
and nyg/nb < lO-^-lO"'*. 



C. Yqq + Yqq (T < Tt.a) 



Q 



Yqq + Yqq O {YYq) + qqq, 
: EsiiYYq)] -mN + 2Am{Yqq). 



(43) 



and large more likely than that of the first reaction [Eq. 
([55])]. and less likely than that of the second [Eq. P^ ]. 

If we refer to values of charmed baryon case, the Q- 
value is negative: The mass of A+ = udc is 2284 MeV. 
The mass of S+ = ucc is 3519 MeV. The mass of 
charm quark is 1.27 GeV. This mass difference leads to 
EB[{YYq)] + Am{Yqq) = -0.979 GeV + 1.01 GeV = 
—0.03 GeV [2^. Adopting this and the assumption, 
AmiYqq) = 0.7, we derive Q ^ -0.27 GeV. 

If the previous two reactions [Eqs. (|38p and (|^^ ] are 
not operative in producing the {YYq) and this reaction 
[Eq. (|43p ] is operative for some reason, final abundances 
of (YYq) or {YYY) and Yq would be roughly of the order 
of Yq abundance at the transition temperature or less, 
i.e., 3n(^YYY)/nb < lO-^-lO"'' or 2n(^YYq)/nb < 10"^- 
lO"'', and nyg/rib < 10"^-10-'*. 



In Cases 1 and 2, the Q- value is given by Q ~ 
EB\{YYq)]~mN + 1.4 GeV. The Q-value is then positive 



VI. CONCLUSIONS AND DISCUSSION 

We have investigated effects of decay properties of res- 
onances made of two long-lived massive colored particles 
(MCPs), i.e., (YY), and binding energies or energy levels 
of exotic strongly interacting massive particles (SIMPs) 
including one MCP in them, on the annihilation of MCPs 
at color confinement temperature Tc ^ 200 MeV. Mag- 
nitudes of cross sections of two-SIMP collision for final 
states including a {YY) is expressed by a parameter f^. 
In addition, binding energies of SIMPs are taken as other 
parameters. We studied effects of the parameters on 
the YY annihilation, and resulting relic abundances of 
SIMPs by calculating with different parameter sets. 

It is assumed that the {YY) bound state formation 
from initial states of two SIMPs dominantly proceeds 
through resonances (or barely-bound states) of {YY). 
The cross section for bound state formation is then for- 
mulated utilizing a well-known expression for resonant 
cross section. We suppose that bound states or reso- 
nances of MCPs decay into more bound states mainly 
through the dipole photon emission rather than hadron 
emissions, based upon an ansatz from particle data on 
heavy quarkonia |28| . In addition, we roughly estimate 
fates of {YY) states from a comparison of the rate for cas- 
cading down energy levels by the spontaneous emission 
and collisions with background photons to that for being 
destroyed by collisions. We then make a toy model of 
the effective cross section for YY annihilation by two- 
SIMP collisions which accounts for only the reaction 
which ends with the YY annihilation without dissociated 
into two SIMPs by background photons. In this model, 
we roughly take into account the transition of SIMP, i.e, 
Yq — >■ Yqq, where g is a light quark. 

Evolutions of SIMP abundances are calculated as a 
function of the temperature for four cases. Different 
parameter sets for three initial abundances and three 
masses of MCP are considered in the calculation. In 
Figs.[3HSlwe show results of abundance evolution at tem- 
perature T < Tc for some cases. Final abundances for 
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all cases are listed in Table IIVI We show that larger 
values leads to smaller final abundances of SIMP, while 
smaller mass gains of SIMPs lead to larger final abun- 
dances. We then suggest that informations on exotic 
heavy hadrons such as their energy levels and partial de- 
cay widths are necessary to estimate relic abundances of 
SIMPs precisely. The calculated relic abundances (Table 
HV)) for respective SIMP species are 2 x lO'^-S x 10"^ 
times that of baryon depending on parameters. This 
large uncertainty stems from uncertainties in level struc- 
tures and decay widths of SIMPs, and the mass and the 
initial abundance at color confinement of MCP. We fi- 
nally comment on possible scenarios in a case that long- 
lived heavy hadrons composed of two Y particles exist. 

Wc note that the relatively small final abundances of 
SIMPs with a huge uncertainty can be important since 
even such small amounts of SIMPs possibly leave observa- 
tional signatures on primordial light element abundances 
through bindings of SIMPs to nuclei and nuclear reac- 
tions of the bound states [H, [s^ . If a long-lived SIMP 
interacts with nucleon by the strength of the potential 
between two nucleons, the binding of the SIMP to a nu- 
cleon triggers formations of heavy nuclei up to carbon 
which are bound to the SIMP. The decay lifetime of the 
SIMP, therefore, should be tx ^ 200 s to avoid heavy el- 
ement synthesis conflicting with observations of old stars 

If the potential between a long-lived SIMP, i.e., X, 
and a nucleon is suitably weaker than that between 
two nucleons, primordial abundances of ^Be and ^Li 
can reduce by the reactions X-|-''Be^'^He-|-*IIex and 
X-|-''Li— >■ t+'^Hex, respectively [13]. The discrepancy 
in ^Li abundance between predictions in standard big 
bang nucleosynthesis model and observations of metal- 
poor stars |35l - [4]| can be caused through the exotic nu- 
clear reactions enabled by the existence of long-lived relic 
SIMPs. The abundance of SIMPs necessary for a reduc- 



tion of ^Li abundance down to the observed level in an 
example case (figures 5 and 7 of Ref . [s^l ) was ~ 1.7x10"^ 
times as large as the baryon abundance. That relatively 
high SIMP abundance is predicted especially in Case 3 in 
this study. The present study indicates that a high relic 
abundance tends to result from small branching ratios 
of (YY) states for decays to lower (YY) bound states, 
a large initial abundance of MCP, and large MCP mass. 
Although we did not treat a case that the colored parti- 
cle Y is electrically neutral, its relic abundance may be 
larger than in charged cases as we have mentioned at the 
end of the Sec. IIIII The neutral Y may, therefore, also be 
a viable candidate for SIMP needed for a ^Li reduction. 

There are interesting possible predictions to be tested 
in future astronomical observations adding to that on ^Li 
[3^. One is ^Be and/or B abundances larger than the 
prediction of standard big bang nucleosynthesis. Another 
is a high isotopic ratio ^°B/^^B different from predictions 
of other models for boron production such as the cosmic 
ray nucleosynthesis (^°B/^^B^ 0.4 [42 - |44| ') and the su- 
pernova neutrino process (^°B/^^B^ 1 [4^, E^). 

If a long-lived SIMP affected the nucleosynthesis, pri- 
mordial light element abundances would have a unique 
information on the quark-hadron phase transition in the 
early Universe. The relic abundance of SIMP signifi- 
cantly depends on hadronic characteristics through an- 
nihilation rates as shown in this study. In view of this, 
more elaborated studies on the annihilation process are 
necessary to derive an information on the phase transi- 
tion in future. 
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